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Abstract - We study semilinear Maxwell-Landau-Lifshitz systems in one space dimension. For 
highly oscillatory and prepared initial data, we construct WKB approximate solutions over long 
times 0(l/e). The leading terms of the WKB solutions solve cubic Schrodinger equations. We show 
that the nonlinear normal form method of Joly, Metivier and Rauch [4] applies to this context. 
This implies that the Schrodinger approximation stays close to the exact solution of Maxwell- 
Landau-Lifshitz over its existence time. In the context of Maxwell-Landau-Lifshitz, this extends 
the analysis of Colin and Lannes [2] from times 0(\ In e|) up to 0(l/e). 
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1 Introduction 

The state of a ferromagnet is described by the magnetization vector M. The evolution of M with 
time t without damping is governed by the Landau-Lifshitz equation 

d t M = —M x H, (1.1) 

where H is the magnetic field. The electromagnetic field variables obey Maxwell equations 

f d t E - V x H = 0, , N 

I ' 1.2 

| d t H + V x E = -dtM. 

The Maxwell-Landau-Lifshitz system that we study is (jl.ip ~ p.2p . The physical constants have 
been set equal to one. The spatial domain is R 3 . Observe that the Maxwell-Landau-Lifshitz 
system admits a family of constant solutions: 

{E,H,M) a = (0,aM ,M o ), 

for any a > and any Mq £ R 3 . We are interested in the solutions that can be written as small, 
slowly variable perturbations of such constant solutions. The perturbations are measured in terms 
of an arbitrarily small positive parameter e. We postulate the following form for the perturbations: 

E(t, x) = sE(et, ex), H(t, x) = aM + sH(et, ex), M(t, x) = M + eM(et, ex). 

If (E, H, M) is a solution of (fL~T j) -(jL~2 ]) . then (E(t, x),H(t, x),M(t, x)) satisfies the following system 

' d t E -V x H = 0, 

< d t H + V x E = -dtM, (1.3) 
k d t M = -e^M x H + e^aM x M - M x H. 



3_iH9_d3 



Introduce the perturbation unknown u = (E, H,a? M) : [0, oo) x M 3 — > 
the system p.3p can be written as a symmetric quadratic hyperbolic system in 



x 



u: 



d t u + A(d x )u + -L u = B(u,u), (1.4) 

e 



where A and Lo are defined by 

/ -d x x 0\ /0 

A(d x )=id x x , L = -M x a-M-x | . (1.5) 
V 0/ \0 a-M x ~aM Xj 

The bilinear form B(-, ■) is defined for any u — (u 1 , u 2 , u 3 ), v — (v 1 , v 2 , v 3 ) £ I 3 x R 3 x K 3 by 

if ° \ 

B(u, v) = - I a a ( u 3 x w 2 + w 3 x u 2 ) J . (1.6) 

\ — (u 3 x v 2 + v 3 x u 2 ) ) 

It is easy to check that A(d x ) — Ajdj with Aj symmetric and that Lq is skew-symmetric. 

We set a = 1 and M = (1, 0, 0). 
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In [7] , Leblond did formal computations for Maxwell-Landau-Lifshitz system, and derived a Davey- 
Stewartson system in the high-frequency limit, for diffractive times. Our goal is to justify the 
Davey-Stewartson approximation over long time interval of order 0(1/ e). For technical reasons 
that are made clear in Section 1.1.3 and have to do with the issue of localizing the resonances, so 
we restrict our problem to one spatial dimension, where the Davey-Stewartson system degenerates 
into cubic Schrodinger equation. In one spatial dimension, the equation (jl.4l) becomes 

dtV + A(e x )d y v+^ = B(v,v), (1.7) 

where the space variable is y £ R 1 , and the vector e\ € K 3 is 

ei = (1,0,0)*. (1.8) 
We consider highly oscillatory initial data of the form 

v(0,y) = a( y y k y/ £ +^y)e-* k y^ +ea 1 (y,ky/e) + e 2 a 2 (y,ky/e), (1.9) 
where ai(y,9) and a 2 (y,9) are real- valued and 27r-periodic in 9. 

The system (|1.7p is semilinear symmetric hyperbolic. For fixed e > 0, the local existence, unique- 
ness and regularity of local-in-time solutions to (11.7I) - (|1.9[) in smooth Sobolev spaces are classical. 



Our goal here is to describe the solutions in the high-frequency limit e — > over long times 0(l/e). 
There are three main issues. 

(1) . First, because of the fast oscillations, the Sobolev norm of the initial datum is unbounded as 
£ — > 0: for any s > 0, 

|«(0)| H . = \a(y)e^ e + ...\ H . = 0( £ - s ) -> oo. 

By the classical theory of Cauchy problem of symmetric hyperbolic system, the existence time 
t*(e) shrinks to as e — ► 0. 

We can go around this issue and show that the natural existence time is 0(1) by introducing 
profiles as in [3] and [4]. This means we consider a representation of the solutions by maps 
depending periodically on the fast variable 9 € T = R/27rZ: 

v(t,y) = V(t,y,8)\ B _ ( - ut + kv) t (1.10) 

where lo € M. is chosen such that the couple (w, k) satisfies dispersion relation det (— iuj + A(e\)ik + 
Lq) = 0. Then to solve the Cauchy problem (|1.7|) - (|1.9|) . it is sufficient to solve the following Cauchy 
problem in V: 

d t V + A(e l )d v V + ^{-ujd e + A( ei )kd e + L }V = B(V, V), 

1/(0, y, 9) = a(y)e ie + ^)e~ td + ea x (y, 9) + e 2 a 2 (y, 9). 

This is a symmetric hyperbolic system, for which the classical theory gives local-in-time well- 
posedness, uniformly in e. Indeed, the data are now uniformly bounded with respect to e in 
Sobolev spaces with non-negative indices, the L 2 estimate is uniform in e in spite of the large 
1/e prefactor because the operator — ujdg + A(e\)kde + L is skew-adjoint, and the commutator 
estimates are trivial because the operator — uide + A(e\)kdg + Lq has constant coefficients. 

(2) . Second, the initial datum (jl.9l) has a large amplitude 0(1). In this regime, the maximal 
existence time a priori is 0(1), not 0(1/ e). For (|1.11|) . the classical H S (M. x T) energy estimate 
for semilinear symmetric hyperbolic operators is 

V(t)\ H ° < \V(0)\ H s+C f \V(t')\ H *dt', 
Jo 
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where the constant C depends on |V|l°°, implying, by Gronwall's lemma, the bound 

\V(t)\ B . < \V(0)\ H .e Ct . 

Over long times 0(1/ e), the upper bound diverges to infinity as e — > 0. We will overcome this 
problem by normal form reduction and a rescalling method, following Joly, Metivier and Rauch 
in [4] and Texier in [9]. This is in contrast to the situation in which the initial data have a small 
amplitude 0(e), where the regime is said to be weakly nonlinear and the natural existence time is 
0(l/e). 

(3). Third, we are confronted with the issue of deriving uniform bounds for the corrector terms of 
the ansatz. If we denote V = V Q + eV 1 + ... the WKB expansion, we want indeed to guarantee 
that eV 1 is o(l) in the time interval under consideration, otherwise the significance of the WKB 
expansion as an asymptotic expansion breaks down. However there typically arises an equation of 
the form (dt + pd y )a = b, where a is a component of V 1 and b depends on the leading term V . 

For the trivial datum a(0, y) = 0, the solution is a(t, y) = b(t — s,y — ps)ds, giving typically the 

Jo 

estimate ||a(t)||jj» < i||6||,L°c(ifs) and nothing better. Over long time 0(1/ e), the upper bound is 
of order 0(l/e) implying that eV 1 could well be 0(1). This phenomenon is called secular growth 
and was extensively studied by Lannes in [BJ. 

For Maxwcll-Landau-Lifshitz, we overcome the issues (2) and (3), and show existence up to times 
of order 0(l/e). We also show that the solution can be well approximated by WKB solutions, the 
leading terms of which satisfy cubic Schrodinger equations. 

In the rest of this introduction, we describe in greater detail the issues associated with (2)-(3) 
above in upcoming Section 1.1, before turning to the statement of our results (in Section 1.2). 

1.1 Context 

1.1.1 Secular growth 

In [2] , Colin and Lannes study systems of the form (|1.4p with highly oscillatory initial data 

u(0, x) = a(x)e fk - x ' e + ^x) e - ltx/e + ea x (x, —) + e 2 a 2 (x, — ) + (1.12) 

e e 

corresponding to (jl.9p in several space dimensions. 

They first show, under a weak transparency assumption (see Assumption 2.1 in 4 ) ensuring that 
a WKB cascade exists that an a three-scale approximate profile of the form 

2 

V a (t,x,0) = 1 £ / E j V j (T,t,x,0)\ T = et , V 3 (r,t,x,6) = J2^ p6 V jp (r,t,x) 

j=0 p 

can be constructed by three-scale WKB expansion over long times of order 0(1/ e). The leading 
terms of the approximate solution solve Davey-Stewartson systems. Under the additional assump- 
tion that the Davey-Stewartson systems are well-posed, this proves existence of WKB solutions 
over time 0(1/ e). As mentioned in Section 1.1, Colin and Lannes have to deal with the issue of 
secular growth. There arises in their analysis an equation of the form 

(d t + u'(fc) • d^irtVn = 27riB(V i, (^10)*), (1-13) 
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where ttq is a spectral projector (see (|2.11|Q and (ttoVio)* is from the following decomposition 



TToVio = (7T Vio)* + (7T Vio) 

where the part (7ToVio) satisfies the transport equations (<9 t +u/(fc) -d x ) (710 Vio) = and the residual 
(tto^io)* satisfies the following homogeneous hyperbolic equation 

(d t + n A(d x )n )(ir V 10 )* = 0. (1.14) 

The couple (V01, (""0^10)) is the solution of a Davey-Stewartson system and is uniformly bounded 
in a Sobolev space over time interval [0,T] r x R t , where the vector (7r Vio) is determined by the 
leading term Vq±. 

Given any initial datum to (7r Vio)*, the equation (|1.14[) implies a global solution (7r Vio)*. Back 
to (|1.13|) . one may only generally have an estimate ~k\V\\ = 0(t), unbounded as t — > 00. Assuming 
that the initial datum a e T, s = {u e H s+n+ / 2 , (1 + 1 • \) s u e for some s > 3 + n/2, Colin and 
Lannes prove a better secular growth estimate in O(lni) for ttxVu (this is Proposition 2.3 in [2]). 
Then, assuming that the initial perturbation (v — v a )(0) is 0(e 2 ), Colin and Lannes prove stability 
of the WKB solution in times 0( \ lne|) (this is Proposition 3.1 in [5]). 

Here, we consider the problem in one space dimension, with initial data in classical Sobolev spaces. 
We allow initial perturbations of size 0(e), and give two results (see Theorem ll.3l below for a precise 
statement), the proofs of which rely heavily on the specific structure of the Maxwell-Landau-Lifshitz 
system. 

First, assuming only the standard polarization condition for the leading amplitude a in the initial 
datum (|1.9j) . that is a(y) G ker(— ico + A(ei)ik + Lq), we give a stability result in intermediate 
times. This is statement ii in Theorem 1.4 below. We show that a WKB approximate solution 
can be constructed, of which the leading term solves cubic Schrodinger equation. We choose 
(ttoVio)*(0) = to eliminate the secular growth of it\V\\. Indeed, with this null initial datum, the 
equation (j!.14|) implies that (7ToVio)* is identically zero. Then the source term in (|1.13j) is null, 
hence no secular growth for "K\V\\. The initial error is of order 0(e): 

\(v-v«)(0)\ L ~=O(el (1.15) 

where the functions v and v a are respectively the exact solution and the approximate solution. We 
then prove the error estimate in time 0(e a \ lne|/e), for any < a < 1: 

\( v - v a)( t )\ =0(e a \lne\). 

Second, for the above WKB solution, we give a stability result in diffractive times 0(l/e), corre- 
sponding to statement iii in Theorem 1.4 below. We assume additionally that the correctors are 
prepared as follows 

\Il s (v - v a )(0)\ L ~ =0{e 2 ), (1.16) 

where the eigen-projectors ±1 and II S are defined later in (jl.29[) . In the context of our WKB 
analysis in Section 2, the condition (j!.16[) is equivalent to the following condition (see also (I2.51[) ): 

ax(y,9) = (1 - Tro)Fio(O) + <^io>(0) + a(y,9), (1.17) 

where a\ is the first order corrector of the initial datum (|1.9[) . and a is a vector function that 
satisfies IT s a = 0. Under the preparation condition (|1.16|) . using the specific structure of Maxwell- 
Landau-Lifshitz (see the upcoming Section 1.1.2), we then prove that, up to times 0(l/e), the 
approximate solution is stable in the following sense 

sup \n (v-v a )\=O(e), sup \U s (v-v a )\ = 0(e 2 ), 

te[0,T/e], V eR *6[0,T/e],yeR 
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meaning that the initial errors (|1.15l) and (I1.16[) are essentially preserved by the dynamics of the 
Maxwell-Landau-Lifshitz system in the dim-active regime. 

The condition (|1.16[) is weaker than the corresponding condition (24) in 2 . Indeed, we only give 
constraints on the component of the first corrector that belongs to the image of II S . Unlike the 
condition (24) in 2 , Vii(O) does not appear in (jl.l6[) . because we have the equality n s yn = 
(see (gUl ). 

1.1.2 Maxwell-Bloch structure 

In PP, Colin (following Joly, Metivier and Rauch [3]) studies systems of the form 

' {d t + Bx(d x ) + -)u = f(u,m), 

p (1.18) 
(d t + B 2 (d x ) + — )m = g(u, u). 

This operator has a block-diagonal structure which is very particular. This is the one which is 
used in [3]. Its characteristic variety is the union of the characteristic varieties of d t + Bi(d x ) + E 
and dt + B2(d x ) + F. Assume the spectral decomposition 

3 3 

where Xj and fa are eigenvalues, nj and Ll| are eigenprojectors. 

Colin makes a strong transparency hypothesis as in Joly, Metivier and Rauch [4] (see Assumption 
2.10 in g]), saying that 

Assumption 1.1. There exists a constant C > such that for all £i,£2>?7 € K d which satisfy 
rj = £i +^2,. f or an U eigenvalues Hj{rj), Aj/(^i), A_y"(^2) respectively of B^{rj) + F/i, i?i(Ci) + E/i 
and Bi(^2) + E/i and for all u, v G C", one has 

|n?(»?)ff(n},(^)u,n)„(&)t;)| < c\\ j ,(t 1 ) + \ jll fa)- »M\\u\\v\. (1.19) 

We say a system has Maxwell-Bloch structure, if it is of the form (|1.18[) and satisfies Assumption 

o 

For systems with Maxwell-Bloch structure, Colin shows that a WKB expansion can be performed 
on (|1.18l) . where the leading terms of the three-scale WKB solution satisfy an elliptic Davey- 
Stewartson system which is locally well-posed, implying that an approximate solution (u a ,m a ) 
can be constructed over times 0(1/ e) (see Proposition 4.1 in [1]). Assuming additionally that the 
initial correctors are identically zero, that is, (u,m)(0) = (u a ,m a )(0) (which is strictly stronger 
than (|1.16|l ). Colin showed a convergence result over times 0(1/ e). The Maxwell-Bloch structure 
gives a control on the interactions of resonant waves (see Section 4.2 in [T] and Section 3.2 in this 
paper). In this respect, it is crucial in the derivation of error estimates in time 0(1/ e). 

Here, we check that the Maxwell-Landau-Lifshitz system in one space dimension (jl.7p is of Maxwell- 
Bloch type. 

1.1.3 Resonances 

In order to check Assumption 1.1 on (|1.4j) . we need to describe the resonances. For any £ = 
(£i> £2> £3)1 w e have the spectral decomposition 
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3 9 
m— 1 j — 1 

The resonance set is defined as the domain in frequency space in which the upper bound in 
vanishes, that is the set of frequencies (£, 77) such that $j 1 J2 J3 (£,77) =0, where the phases functions 
are defined by 

*iiJaJ»K.»7) = A 3i(£ + ^ - A 3 2 (0 - Aj.fa)- 

For (jl.4l) , the eigen-polynomial is 

dct ( - A + J2 Mi + Lo/i) = A 3 [A 6 - 2(2 + |£| 2 )A 4 + (|£|> + |£| 2 ) - 2£ 2 )A 2 - |£| 2 (2|£| 2 - &)] . 

3=1 

It is of degree nine with three parameters £1, £2 and £3. The eigenvalues Xj do not symmetrically 
depend on £1, £2 and £3. This property makes the resonance equation &j lt j 2 ,j 3 (£, Tj) = too 
complicated to solve analytically. For this technical reason, we consider the one space dimensional 
system (jl.7|) . in which the dispersion relation reduces to a simpler formula (see below (|3.3l0 . 

1.2 Description of the results 

If a C 1 profile V(t, y, 0) solves the Cauchy problem 

8 t V + ~A{e x )d v V + - £ {-Lod e + A( ei )kd e + L }V = B(V, V), 

V(0, y, 0) = a(y)e te + ^)e- ie + eai{y, 0) + e 2 a 2 (y, 0) 

on the time interval [0,T/e], then the function v(t,y) defined as (| 1 . 10[) solves the Cauchy problem 
([L7 |) -(fL9 )l on the time interval [0,T/e]. 



Throughout the paper, we always make the following assumptions: 

We assume that the phase (w,fc) satisfies w ^ and the dispersion relation det L(i(ui,k)) = 
with L(i(uj, k)) :— {—iuj + A(e{)ik + Lq) = 0. This means that we choose a time frequency that is 
compatible with the initial space frequency. We assume that the leading term of the initial datum 
satisfies the polarization condition: a(y) £ kerL(z(w, k)). We also assume the regularity conditions 
a(y) £ H s , ai(y,0) £ H l (T e , H^ 1 ^)) and a 2 (y,0) £ H^Tg, iP- 2 (K*)), where s > 2 + 1/2. 

In Section 2, we construct an approximate profile for (|1.20p by WKB expansion. This construction 
relies on the structural condition (|2.9I) . which implies the weak transparency condition of Joly, 
Metivier and Rauch. 

Proposition 1.2. (Existence of Approximate Profile) 

Under the above assumptions, there exists V a £ L°°([0, T/e] t , ^(Tg, iJ^ 2 (R^))) for some positive 
T > independent ofe, such thatV a solves the following Cauchy problem on time interval [0,T/e] t : 

d t V a + A( ei )d v V a + - £ {-Lod B + A( ei )kd e + L }V a = B{V a , V a ) + e 2 R 
V a (0, y, 0) = V(0, y, 0) + eb(y, 0) + e 2 b x (y, 0) 

for some R(t,y,9) £ L°°([0, T/e] t , H x (T e , ff- 2 (!j))), {b{y,9)MvM & H\T e , H s - 2 (Rl)). 
Precisely, the function b is given in (|2.47p and has the property stated in Lemma 12. U the function 
R is given by (|2.49[) and has the property stated in Lemma 12.21 Moreover, the leading term of the 
WKB approximate profile V a solves a cubic Schrddinger equation. 
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We remark that it does not seem possible to construct a more precise WKB solution. Indeed, for 
higher order terms, there arises an equation of the form 

(-p + TT A(e 1 )7r )d y (V 20 ) =4Re7roS(Fo,- 1 ,F 2 i) + 2Re7r B(F 1) _ 1 ,F u ) -d T (V w ), (1.22) 

where the right hand side is not of the form d y lZ, then equation (ll.22j) is not well posed in classical 
Sobolev spaces. This can be contrasted with corresponding equations for lower order terms (see 
(|2.28p ). where the source terms are indeed of the form d y lZ. 

For V a given by Proposition HH defining v a := V a (t,x, k y~ ut \ thcn v a £ L°°(\0,T/e}t x Rh) 

S y 

and satisfies the following Cauchy problem on time interval [0,T/e]: 

' Lnt" = U(r'\ ) - :-//. 

(1.23) 



(d t + A(e l )d y )v a + -L v a = B(v a ,v a ) + s 2 R, 



v a (0, y) = v(0, y) + eb(y, ^) + s 2 b x (y, 

s s 

We then justify the WKB expansion by giving an upper bound for \v — v a \ that is valid over the 
existence time for v a : 

Theorem 1.3. Under the assumptions stated just above Provosition 11.21 there exists T* > and 
< Eq < 1, such that for all T < T* , for some constant C — C(s,T) > depending on s through 
\\ a \\fi° + IKHiP( Xj , ;i? s-i) + \\a 2 \\H 2 (T e ,H>>-2)> we have that for all < e < e : 

i. Over the time interval [0,T/s], the Cauchy problem (|1.7[) - (|1.9[) admits a unique solution v of 
the form v(t, x) = V(t, x, ky ~ ut \ with yn Xj Q) e L°°([0, T/e] t , H 1 ^, H s - 2 (Rl))). 

ii. We have the following error estimates 

l|Ho(« ~ v a )\\L-ao^]x^) < C(e + T), ||n a (w - v a )\\ L ^ {[0 T ]xRl) < Cs, (1.24) 

iii. Assuming additionally that ai satisfies (|1.17|) . we have the better error estimates 

\\IL (v - v a )\\ Loa([0 z ]xRl) < Cs, \\U,(v - v a )\\ L ~ {[0 z ]xK i) < Cs 2 . (1.25) 

We remark that for times just short of the diffractive time scale, that is T = s a \ lne|,0 < a < 1, 
the error estimate (| 1 . 24[) gives an upper bound that is o(l) in the limit e — >• : 

ll^-« a )llL~([o, e - 1+ »|in £ |]x R1 ) =0(£ Q |lne|). (1.26) 



The proof of Theorem (|1.3[) is decomposed into two steps. First in Section 3, we show that the 
system (|1.20j) has the Maxwell-Bloch structure, then by a nonlinear change of variable introduced 
by Joly, Metivier and Rauch in [3], we prove that the Cauchy problem (|1.20l) admits a unique 
solution over times of order 0(1/ e). Then in Section 4, we prove the error estimates (| 1 . 24[) and 

We conclude this introduction by setting up some notations: 

For all (A, £) £ R 1+1 , we introduce the matrices L(i(X, £)) = i\ + A(e\)i^+ Lq and the characteristic 
variety Char£(i(A, £)) = {(A,0|det L(i(\£)) =0}. 

For (3 = (u),k), denote L(j3de) = —ujdg + A(ei)kdg + L and L p — L(ip(i) — —ipuj + ipkA(ei) + L . 
Denote by n p the orthogonal projector onto kev L p and L~ x the inverse of L Vl or, if L p is singular, 
its partial inverse which is defined as 

K p L p 1 = L p 1 ir p = 0, L p L p 1 = L p 1 L p = Id — ir p . 
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For all ^Gl'jWe have the spectral decomposition 



3=1 

We denote the total eigenprojectors 

6 9 

n (0 = X)ni(0, n.(0 = Yl u r (0- (1.28) 

3=1 3'=7 

By direct calculation, for all £ G R, we have Xr(^) — A§(£) = Ag(£) =0 and 

Do(0 = diag{0, 1, 1, 0, 1, 1, 0, 1, 1}, n s (0 = diag{l, 0, 0, 1, 0, 0, 1, 0, 0}. (1.29) 

The projector IIo(£) and II S (£) are actually constant matrices, independent of the variable £. 

In the following, most functions and operators depend on e. In this context, we typically use 
bounded to mean bounded uniformly in e £ (0,1). Later on, we restrict the range of e to < e < Eg, 
for some Sq small enough. 

2 WKB expansion and approximate profile 

The aim of this section is to construct an approximate profile of 1(1.200 through WKB expansion. 
In the nonlinear regime of our interest, the limit equations are cubic Schrodinger equations. 

2.1 The WKB Expansion 

We look for an approximate profile V a of the form 

2 

V a (t,y,6) = J2e j V j (et,t,y,6), Vj(r,t,y,e) = ^^(r,^). (2.1) 

3=o pez 

We plug ((2.1[) into (II .20)) and cancel all the terms of orders 0(e J ), j = —2,-1,0. 

Equations for the terms in 0(e" 2 ): L(/3de)V a = (-cudg + A(ei)kdg + L )V = 0. By l(2T]) . it 
amounts to solving the following equations for all p S Z: 

LpVop = (-ipuj + A{ ei )ipk + L )V 0p = 0. (2.2) 

We study for each p E Z. 



p=0 If p = 0, the equation (|2.2[) becomes LqVqo — 0. Since there is no mean mode in the leading 
term of the initial datum (11.200 . we choose naturally Vbo = 0. Actually, assuming the leading mean 
mode Voo to be null is crucial to obtain a Davey-Stewartson approximation, as in [I] and [2] (here 
the Davey-Stewartson approximation degenerates into a nonlinear Schrodinger approximation). 



p=l If p = 1, the equation ((2.21) becomes LiVqi = 0. By direct calculation, the dispersion 



relation det L\ = is 



^ - 5> a = < 2 - 5> a - ^ 
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We choose lu ^ and k ^ that satisfy the dispersion relation (I2.3[) . We denote 

S = -^=#l, T .I-4, ,2.4) 
then the dispersion relation (|2.3[) is equivalent to 

2 W + 2(5 2 
K = — • U) , = ±1. 

Through direct calculation, we find that ker L\ is a one-dimensional vector space with generator 

Wo = ft , fio = W • (2.5) 

V-7O0/ V 1 / 

Then the orthogonal projector tt± onto ker L\ has the expression: 

7T!V= \ ' °V , for all FeM 9 . (2.6) 

I vvo r 

where (-|-) denotes the inner product in the complex vector space C 9 . Then the solution Vbi has 
the following form: for some scalar function g (which is determined later by (|2.40p and (12.411) ): 

V01 = W09 (2.7) 



|p| > 2 The equation (|2.2p becomes L p Vo p = 0. For nonzero (uj,k), det L p = -^=4> p 2 w 2 (l — 

jjz) 2 = (2 — ^-) 2 . By our choice of (w, fc) that satisfies (12.31) . the matrix L p is invertible if \p\ > 2. 
Then the solutions are trivial Vo p = 0, for all \p\ > 2. 



p < — 1 By reality of the datum, the only natural choice is Vj tP = V Jt - p . 



Equations for the terms in 0(e 1 ). For all p: 



(d t + A(e 1 )d v )V 0p + L p V lp = J2 B ( V 0p^V 0p2 ). (2.8) 

Pl+P2=P 



p—0 If p = 0, the equation (|2 . 8[) reduce to 

L V 10 = 2B(V 0U V ,-i). 

We calculate 

B(Voi,Vo,-i) = ggB(W ,W ) = y ( -7^0 x_po - Tpo x fi | = 0. 

\ 7^0 x fi + 7^0 x fl 

This means we have B(wi-,%_x-) — 0- The relationship noB(wi-, tt-%-) = 0, which is necessary to 
construct the Davey-Stewartson approximate solution in Colin and Lannes [2], is the weak trans- 
parency condition introduced by Joly, Metivier and Rauch [4]. Here, we even have the following 
stronger result 

B(7r±r,7r± r )=0. (2.9) 
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Then the equation (|2.9p becomes LqViq — 0. We denote ttq the projection onto kerio- Then 
we have Via — ttqVxo. Now we calculate the vector space kerLo and the projector ttq- We let 
Vio = (Eio,i?io,Mio)*, where E w , H w and M w are vectors in R 3 . The equations (|1.5|l and (|2.9|) 
give 

L V W = M x (ff 10 - Mxo) = (fr M - Mio) // (1, 0, 0)*. (2.10) 
Then we obtain a basis of kerL : 

ker L — span{ I I , 

w 

where e\ = (1, 0, 0)*, e 2 = (0, 1, 0)*, e 3 = (0, 0, 1)*. The orthogonal projector 7r onto kerL is 

7T = [ Jr J 2 I (2.11) 

with the blocks 








Jl = i , J 2 = i . (2.12) 



p=l If p = 1, the equation becomes (d t + A(ei)d y )V i + LiVu = J2 Pl + P2 =i b ( v Qpi > v o P2 )- 



From the analysis for the terms in 0(s 2 ), we have that Vbo = and Vo p = 0, for all \p\ > 2, which 
give us J2 Pl + P2 =i B(V 0pi , V 0p2 ) = 0. Then the equation for p = 1 is 

(ft + A( ei )d y )V 01 + LtVu = 0. (2.13) 

Applying m to (|2.13|) . by identity "K\L\ = 0, we obtain 

iti(d t + A{e x )dy)V m = Q, (2.14) 

a transport equation for the leading profile. In the geometric optics approximation, the transport 
occurs at group velocity (See Proposition 2.6 of [5] and Proposition 2.2 in [5] ). Here we compute 
the transport operator -KiA(e{)d y 'K 1 explicitly. By (12.61) . the equation (12. 14)) is equivalent to 

((d t +A(d x ))V ai \W Q )=0. 

By (1231 and (|2~7)) . we have that 

(d t V Q1 \W ) = d t g |Wo| 2 = d t g 2(l + ^ + (I - ^f) 

uj z 

and that 

4k 

(A(e 1 )d y Voi\W ) = —d y g. 
Then the equation (|2.14p becomes the following transport equation in the scalar function g: 

d t g + pd y g = (2.15) 
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with the group velocity 

(k , k . t-t\ — 1 2fc 

1 + — + 1-— 2 -• 2.16 

On the other hand, by applying L^ 1 to (|2.13[) . one has 

{l-TT 1 )V 11 = -L^ 1 ((dt+A{e 1 )dy)V m )=-L^ 1 (A{e 1 )d y V Q1 ). (2.17) 
Explicitly, the solution of equation (j2 . 13[) is Vn — [E\\, H\\, M\\f with 



(2.18) 



En = 


--mo 

u 


+ -(-1- 

a; 


- — 

Ul 


< #n = 








Mu = 

< 


-7/fio + 


2ifc, 

— -U 


kp 
u 



where / is an unknown function that is yet to be determined (see ()2.43|) ). and is associated with 
the component itiVn of Vy\. 



p=2 If p = 2, the equation becomes (d t + A(ei)d y )V 02 + L 2 Vi 2 = J2 Pl + P2=2 B ( v o P i , V 0p2 ). 



Since Vbo = 0, Vq p — 0, for all \p\ > 2 and (|2.9[) , the equation for p — 2 is actually £2^12 = 
S(Voi, Vqi) = 0, which gives trivial solution V12 = because the matrix L 2 is invertible. 



p>3 The equations are (d t + A(ei)d y )V 0p + L p Vi p = J2 Pl+P2 = p B (V 0pi , V 0p2 ), and we directly 
obtain the solutions Vi p = 0. 



p<0 For negative p, we simply take Vj P = Vj t - p . 
Equations for the terms in 0(e°). For all p, 

d T V 0p + (d t + A( ei )d y )V lp + L p V 2p = 2 B(V 0pi ,Vi P2 ). (2.19) 



P l+ P 2 = 



p=0 For p=0, the equation (|2.19|) becomes 



d T V o + (dt + A(e 1 )d y )Vio + L V 2 o=2 £ B(V 0pi , V lp2 ). (2.20) 

Pl+P2=0 

We calculate the right hand side 

]T B(Vo Pl ,Vi P2 ) = B(V 01 ,V 1 ,- 1 ) + B{V ,- 1 ,V 11 ) = 2ReB(Vo,-i,V 11 ). 

P l+ P 2=0 

Together with the choice Vbo = 0, the equation (|2.20p becomes 

(d t + A{ei)d y )V w + L V 20 = 4Refl(V 0l -i ) V n ). (2.21) 
Applying ttq to (|2.21l) . we obtain that 

(fit + 7r i4(ei)^,7ro)7r Vio = 47r ReB(V ,-i, Vn). (2.22) 
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By ([H]), (J23]), dH7||, (I^TTj) and (|^T5j) . we calculate the source term of (l2~22l : 

47roHaB(7o,_i,7u) = 

with 




A = 

By {T3J}, ([2TT1) and we have that 

7r A(ei)7ro 



fc 



4fc<5 

CO 



— (l-^)W 2 



(-eix)J! (-eix)J 2 
Ji(eix) 

J 2 (eix) 



(2.23) 



(2.24) 



(2.25) 



We denote V w = (_E 10) #io, Mio), by (ET2"3"1) and (pl2"5l) . the equation (pl2"2"j) becomes 

' a t ^i - (eix)Ji^iJio - (eix)J 2 d y M w = 0, 

< a t ffi + J x {e x x)d y E 1Q = A , (2.26) 

t d t M 1Q + J 2 (e 1 x)d y E 10 = -A Q . 

We now introduce a decomposition inspired from the averaging method of Lannes |5 5 : for any vector 
function V(t, y), we use the notation (V) and V* to denote the two parts of the decomposition 
V = (V) + V* , where the part (V) satisfies the transport equation (<9 4 + pd y )(V) — with the 
same group velocity p as the leading term Vol. To solve (12.26[) . it is sufficient to solve the following 
two systems: 

- pd y {E w ) - (eiX)Jid y (H w ) - (e 1 x)J 2 d y {M w ) = 0, 

- pd y {H w ) + J 1 (e 1 x)d y (E 10 ) = A , (2.27) 

- pd y (M 10 } + J 2 {eix)dy(E 10 } = -A Q 
and 

8 t E* w - {e x x)J x d y Ht Q - { ei x)J 2 d y M* w = 0, 
d t H$ Q + J x {e x x)d y E* w =Q, 
d t M* + J 2 (e x x)d v E* 10 = 0. 
By taking — pd y to the second and third equations of (|2.27l) . we obtain that 

f p 2 d 2 y (H w ) + J 1 {e x x)[{e x x)J x d 2 y {H w ) + (e 1 x)J 2 d 2 y (M w )] = -pd y A Q , 
\ p 2 d 2 y {M w ) + J 2 {e 1 x)[{e 1 x)J 1 d 2 y {H w ) + (e x x) J 2 d 2 y {M w )] = pd y A . 

Since Vio = {E X q, i?io, M xo ) G ker L , by (|2.10[) , we may suppose that 



(2.28) 



(2.29) 



Hi 



Mi 




(2.30) 



Then by (L8\i and (|2~T2)l . the equation ((2391) becomes 



-(1 - —)d v \g\ , 



{p 2 -^)d 2 y (h 2 ) = { P 2 --)d 2 y(hz)=Q, 



(2.31) 
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A solution to (12.311) is 

_4Hp fcp 4H /cp = = Q 

cj w w 

Plugging (|2.30[) and ()2.32|) into the first equation of (|2.27|) . we have that d y (Eio) = which gives 
the trivial solution (Eio) — 0. 

We solve the equation (|2.28|) by choosing the trivial solution = H^ = M* = 0. By doing so, 
we will not see a source term in (12.38)) and there will be no secular growth for -n{V\\. 

With our choice for and (Vio), we have a solution to (|2.22|) : 

r l( ,= ! //,n I = —(1-^) I . , I W- (2.:«) 
M W J 

where the vector ei is given in (jl.8l) . 

Back to by (12~T)) . ((2~Tgf and ([2331 . we obtain that £0^20 = 0, which admits the trivial 

solution V20 = 0. Now equation (|2.21l) is solved. 





p=l For p=l, the equation (|2.19[) becomes 

d T V 01 + (d t + A( ei ))d v V u + L X V 21 = 2B(V ai , V w ). (2.34) 
Applying wi to (|2.34|) gives 

d T V 01 +7ri(9 t +A{e 1 ))dyV n = 2^B{V m ,V w ). (2.35) 

By (j2T7)) . we have that 

n^dt + Aie^dyVu = (d t + tti A( ei )fi^ri Ki^i " ^A{e x )d y L^ (A{e x )d y V 01 ) , 
We already have iriA{ei)d y 'Ki = pd y with group velocity p in (|2.16[) . In diffractive optics ap- 
proximation, we have — niA(ei)d y L^ 1 A(ei)d y — — —ui"(k)dy, which gives a Schrodinger equation. 

(See Proposition 2.6 in [5] or Proposition 4.1 in |3J. Here we compute the Schrodinger operator 
— niA(ei)d y L^ 1 A(ei)d y explicitly. 

By (12.61) . the equation (12.35[) is equivalent to 

(d T Voi\Wo) + ((d t + A(e 1 )d y )V 11 \W ) = 2{B(V 01 ,V w )\Wo). (2.36) 
Since V 01 = gW , we have (d T V 01 \W ) = d T g\W a \ 2 . By ([23]), (f2~15|) and (j2~18l) . we have 

(d t Vu\W ) = \W \ 2 d t f + ^(1 - — )(1 - 2 7 )9 y 2 ff 

and 

(^(exJ^nlWb) = -d y f - -(1 - ^)d> 
For the right hand side of (j2~36) . by (HU), (J2U) and (j2~53")) . 

2(B(V 01 ,V 10 )\W ) = ~ k A{l 1 2 )9\9\ 2 - 

up UJ 
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We denote the real constant 



1 +- + ( 1 -^) 2 ' 



2 oj 2 v cj 

then the equation (|2.35[) is equivalent to 

9r<7 + -[^(1 - 2 7 ) - 1](1 - k Adl9 + (d t + pd y )f = —(1 - - 7 2 ) ff | 5 | 2 . (2.37) 

z/w oj w 

We decompose (|2.37l) into two equations: a transport equation in / related to ttiVh: 

(d t +pd y )f = 0, (2.38) 
and a cubic Schrodinger equation in g related to Voi : 

d r g + iv^g = iv 2 g\g\ 2 , (2.39) 
where the real constants v\ and v 2 are defined as 

* = ~ [^(1 - 27) - 1] (1 - ^ = -*£)(! - 7 2 )- 

a; a; w 

Since the scalar function g satisfies the transport equation (|2.15[) . then g is of the form 

g(r,t,y)= gi {T,z)\ z= y- pU (2.40) 

for some scalar function g\. By (|2.39[) . the function g\{r, z) solves 

d T gi +iv\d 2 z gi = iv 2 gi\gi\ 2 . (2-41) 

By the classical theory of semilinear Schrodinger equations, given any regular initial datum g% (r, z) = 
ao(z) 6 H s , s > 2 + 1/2, the equation (|2.4ip admits a unique local-in-time solution, with existence 
time Tj* > independent of e, depending on the H s norm of ao- One has the following estimate 
for all T <T*: 

sup ||<?i(r,-)||ff» < C(T) < +00. 

0<t<T 

Then there exists a unique solution g E L°°([0, T*[ T xt f , H s ) to (|2~T5|) and (|239| which takes the 
form (pOU|) . 

Back to equation (|2.34[) . one can obtain the solution V21 = (E21, H2I) M21) in terms of the functions 
/, g and an unknown scalar function h: 

i/S k $ k 

E21 = Mlo + -(-1 + —)d v fn + iua&gSlo, 

oj oj oj y 

H 21 = hn , (2.42) 
Mn = -7^0 + ^(-1 + ^)d y fn a + u 4 d 2 y gSl . 

0J A OJ " 

The constants 1713 and 77J4 are real and defined as 
We define h and / by 

h = 0, f = i[™(-l + h£ ) ]- l UAd g . (2.43) 
L w 2 oj 



15 



Because the function g satisfies the transport equation (|2.15l) . the function / defined by (|2.43|) is 
a solution to (|2.38|) . The point is that, with this choice of h and /, there holds Lemma l2~2l below. 
This is the key property that, together with Lemma f2.f [ allows the change of variable (14. f 2)) in 
Section 4.2. By (pT32"]) and (1213)) . we directly have 

H 2X = M n = 0. (2.44) 



p=2 For p = 2, the equation (|2~T^1) becomes L 2 V 2 2 = 2B(V Q1 ,V 11 ). By (j!T7]) and (|2~T51) . we have 



B(Vqi, Vii) = 0, and by the invertibility of L 2 , the unique solution is V 22 = 0. 



p>3 For p > 3, to equation ([235]) : V^ p = 0. 



2.2 The approximate solution and the remainder 

The vector space kerL(i(w, k)) is one-dimensional with generator Wo defined in (12.51) . The initial 
leading amplitude a(y) £ kerL(i(w, k)) and a £ H s with s > 2 + 1/2, so there exists a scalar 
function a (y) £ H s such that a(y) = a (y)Wo. 

Given initial datum g(Q, 0,y) = ao(y), the transport equation f|2 . 15[) and the Schrodinger equation 
f|2 .39[) admit a unique solution which takes the form (|2.40l) over time interval [0, T*) , where T-j* > 
independent of e is the existence time of the cubic Schrodinger equation (|2.40[) . For any T < T* 
g £ L°°([0,T] T x R t ,H s ). We then choose (f,h) as in (|Q5]l . one has / £ L°°([0,T] T x Rt,^" 1 ) 
and the following estimates for any T <T*: 

By (E2J), we have V 01 £ L°°([0,T] T x R tl H s ). By 1(215)1 . we have y n G L°°([0,T] T x Ej,^" 1 ). 
By (123511 . we have V w £ L oo ([0,T] T x R t ,H s ). By (12^21 . we have V 21 £ L°°([0,T] T x R t ,H s ~ 2 ). 

We then let _ 

V°(t, y, 9) := Vqi (et, t, y)e l9 + V 01 (et, t, y)e^ e , 

■ = {V w + Vne 19 + Vue- W ) (et, t, y), (2.45) 
:= (V 2 ie ie + V 2 ie- ie )(et,t,y). 



VHt,y,t 
V 2 (t,y,( 



Define V a (t,y,9) as 



V a = V" +EV 1 +e 2 V 2 , 



then the profile V a satisfies 



1 



8 t V a + A{e 1 )d y V a + -{-Lode + A{ei)kd e + L }V a = B(V a , V a ) + e z R, 
V a (0, y, 0) = V(0, y, 9) + eb(y, 6) + e 2 bi{y, 6) 



(2.46) 



over long time interval [0, T^/e[. The initial perturbations b and b\ have the expressions 

b{y,e) = -ai(y,6) + (Vio + Vue* 6 + V 1 ie- ie )(0,0,y), (2.47) 
bi(y,6) = -a 2 (y,6) + (V 21 e 10 + V 21 e^ e )(0,0,y), (2.48) 
and the remainder R(t, y, 9) is defined as 

R = -2B(V°, V 2 ) - B{V\ V 1 ) - e2B(V 1 ,V 2 ) - e 2 B(V 2 ,V 2 ). (2.49) 
For any T < T-f , we have the following estimates 

R £ L oo ([0,T/e] t ,£r 1 Cir flj fr*- 2 (Ki))) ) (b,bi) £ ff^Ts.fr- 2 ^)). 
At this stage, Proposition If .21 is proved. 

We show more properties of R and b in the following two lemmas, according to the WKB expansion 
in Section 2.f . 
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Lemma 2.1. For initial perturbations b and b\ defined in (I2.47|) . we have the equivalence 

|n s (v(o) - v a {o))\Hip e ,H-»oq)) = o(e 2 ) n s 6 = o. (2.50) 

Moreover, there holds H s b — i/ and onZy z/ i/ie initial corrector a\ in (|1.9p is given by (|2 . 5 1 1) 

Proo/. By (T09")) and ([2"7Tgj) . we have that II s Vii = 0. Then U s b = 0, by (gUTj) , it is necessary and 
sufficient to have 

4fc<5 ko I ® \ 

ai(y,6) = V w (0,0,y) + MV>6) = (1-— ) -ei \a (y)\ 2 + ~ ai (y,6), (2.51) 

where a% satisfies Il s ai = 0. The projector Il s is defined in (|1.29[) . Direct calculation gives the 
equivalence f)2.50[) 

□ 

Lemma 2.2. With f and h given by (|2.43l) . the remainder R defined in (|2.49|) satisfies 

I1 S R = 0(£) in£ oo ([0,T/e] t) fr 1 (T fl ,H- 2 (Ki))) ) for any T < T*. 
Proof. By (|2.45|) . we have 

V 2 ) = B(Vu,F 21 ) + e 2i9 B(V n ,V 21 ) + c.c. 

By (|2.44p . a consequence of (|2.43l) . the components H21 and M21 of V21 satisfy H21 = -M21 = 0. 
Then by the definition of B in ([T6]l . it is easy to obtain V" 2 ) = 0. Also by (|2~45| . we have 

that 

5(V\ V 1 ) = B(V 10 , F 10 ) + B(Vu,Vu) + e W B(V w , Vu) + e 128 B(Vu, Vu) + c.c. 
By (Trnj), (|2T7|) . (|2~Tgl) and (|!03"]) . direct calculation gives 

i'din.r,,! = ii(r u .r u i = o. i^r, ,.",,)= Mi 1 



where 

Ai= *k_ { kp_ l){ j dyg + fdy g ){ilo x n ). 



With / given in (|2.43p . it is easy to check that fd y g+fd y g — 0, and then B(Vu, Vu) = 0. Through 
direct calculation, we have IL S B{V 10 , Vu) = 0. Then H^V 1 , V 1 ) = 0, and by the definition of R 
in (|2.49|) . we have 

n s i?= - £ n s (2 J B(V 1 , V 2 ) + eB{V 2 ,V 2 )). (2.52) 

The lemma is proved. 

□ 

We constructed a WKB solution v a to (11.23|) . We now investigate the question whether the WKB 
solution actually approximates the exact solution over an interval of existence. 
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Figure 1: The characteristic variety. 



3 Maxwell-Bloch structure and long time exis- 
tence. 



In this section, we first show that the system (|1.20p has the Maxwell-Bloch structure, then prove 
a long time existence result by normal form reduction. 

3.1 Spectral decomposition 

As mentioned in Section 1.2, for the symmetric matrix A(ei)^+Lo/i, we have the following spectral 
decomposition 



The characteristic variety (that is, the union of the graphs £ M- Aj(£)) is pictured on figure 1. We 
have that for any £ e E, 1 < j < 6, and 7 < f < 9: 




(3.1) 



Ai(0>2, A 2 (0>0, 0<Aa(O<l. 



Ar-i(0. V(0 = 0. 



(3.2) 



For 1 < j < 6, we have the dispersion relations 



M0 + 2(S ; 



Sj = (-1)*. 



(3.3) 



In order to reveal the Maxwell-Bloch structure, we define 



6 




(3.4) 



3=7 
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where the total projectors IIo and II S are defined in (|1.28[) and the notation V :— d/i. By (|1.29[) . 
the projector IIo and IT are actually constant matrices. By equation in V, the system in 

(Vb,F s )is 

d t V + -AoVo - —V = I1 B(V + V„, V + V s ), 



d 

d t V s - —V s = n s B(V + v s ,v + v s ), 

£ 

where the Fourier multiplier 

6 

A Q = ^2 X ]( £V v + kVe)n 3 (eV y + kV e ). 
j'=i 

By ((Oil and (fQ9| . we have 

n B(U ,U Q ) =U Q B(n s ,U s ) =U S B(U S ,U S ) ^U S B(U ,U S ) = 0. 
Then the equation (|3 . 5|) becomes 

d t V + -A V - — V = 2U Q B(V 0: V s ), 
e e 

d t V s -^V s =U s B(V a ,V ). 

e 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



3.2 Strong transparency and normal form reduction 

The following proposition states that the system p.8j) satisfies a strong transparency hypothesis. 

Proposition 3.1. There exists a constant C > such that for all £, r\ € M, for any eigenvalues 
Aj(£), Aj'(ry), 1 < j, j' < 6 ; and /or all u, v € C 9 , one /las 

|n s B(n, n,v ( v ) v ) \ < c\x 3 (£) + \ r ( v ) \\u\\v\. 

Proof. By direct calculation, for any 1 < j < 6, the eigenprojector IL, has the formula 



(3.9) 
(3.10) 



where 



ft, 



o 

1 



/ I i Tj = 1 - ^2 > ^ 



and the couple (Aj,£) satisfies the dispersion relation Q3.3p . We then compute 





ei 



(3.11) 



(3.12) 



By p. lip , we have 



7j(0 -Ij'iv) 



(1 



(1 



_ A J -(g) + 2J J - _ Aj/ (77) + 28 y 
1 A,(£) + V A*'(»7)+*i' 



{H0 + Sj){h'{v) + ?>j-Y 



(3.13) 
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Recall that Sj and Sj> can only take the value 1 or —1. If 6j> and 5j have the same sign, by (|3.12[) . 
we have H s B(Qj(£), Qj'{rj)) = 0; if 8ji and 8j have the opposite sign, by (|3. 131) . we have 

By (13 . 10[) and direct calculation, the equation (|3.9|) follows. □ 

With p.8p and Proposition ^. 1[ we see that system [3~51 has the Maxwell-Bloch structure (terminol- 
ogy introduced in Section 1.1.2). This allows to use a nonlinear change of variable introduced by 
Joly, Metivier and Rauch 0]. Together with the preparation condition IT 1^(0) = 0(e), this gives 
existence in time 0(1/ e) for the solution to the Cauchy problem (|3.8p . First, changing variable 
V s = eW s gives 

d t V + -AoV - —V = 2eU B(V ,W s ), 

e e 

O 1 

d t W s - —W s = -U S B(V ,V ). 

e e 

Then we introduce the nonlinear change of variable 



(3.15) 



N = W.-J(V ,V ), (3.16) 
where the symmetric bilinear form J has the following form 

j(E u P eip6 >Y, v « eiqe ) =E ■W«p.««)e <(p+9) '> ( 3 - 17 ) 

p£Z q£Z p,q 

for some J pq to be determined below. The equation in N is 

dtN .^L N = dtWs _ ^l Ws _ j (dtVo _ ^L Vlh Vo) - J(V , d t V - ^>v ) 

e e e e ' 

1 ./ i 



n s B(^ , V )-J{--J2 x j( £V v + kVe)Uj{eV y + kV g )V ,V ) 

. 6 

J(V a , --^2 X 3 {eV y + kVg^ieVy + kV e )V ) - 2j(2eU a B(V , W S ),V ) . 



We choose the bilinear operator J to eliminate the singular term of order 0(1/ e) in the above 
equation. We consider the following equation 



n s B(V , V ) = -U(J2Xj(eV y + fcD^n,- (e2>„ + kV e )V ,V ) 

3=1 

6 

- iJ(V , H eV v + kVe)nj(eV y + kV e )V Q ). 



Equivalently, in Fourier, for all (£,rj,p,q,j,j')'- 

n s B(n,(ery + kp)V 0p ,U r (e(( - rf) + kq)V 0q ) 

= -i\j{er] + kp)J pq (U j (er] + kp)V 0p , Ify(e(£ - rj) + kq)V 0q ) 

- i\j (e(£ -v) + kq),J pq (U J (er 1 + kp)V 0p , U f (e(£ - 77) + kq)V 0q ). 
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A solution to the above equation is given by 

^(n^a.n,,^) :=itt^^|j^p> foralla,fceC*. (3.18) 

By PropositionO J is well defined and is bounded from H Sl (Tg ,H s y 2 )xH Sl (T e , H s y 2 ) to H Sl (Tg , H s y 2 ) 
provided s\ > 1/2, s 2 > 1/2. The system in (Vo, N) is now 

d t V + -A Vo - —Vo = 2eI1 B(Vo,N + J(V ,V )), 

£ £ (3.19) 

d t N - = -2eJ(n B(V ,N + J(V , V Q )),V ). 

e 

We then rescale the time and define (Vo, A/")(t, y, 9) — (Vo,N)(r/e,y,9). The system in (Vo,AT) is 

d T V + \AoVo - ^V = 2U B(V ,M + J(Vo, Vo)), 

£ £ (3.20) 

d T N- ^Af= -2J(U B(V ,Af + J(Vo,Vo)),V ). 

The initial datum is 

V (0, y, 9) = n V(0, y, 9) = ll (e ie a(y) + e~ l9 ^y) + e ai (y, 9) + e 2 a 2 (y, 8)), 

1 (3-21) 

Af(Q, y, 9) = W s (0, y, 8) - J(V , V o )(0, y, 8) = -U S V(Q, y, 8) - J(V , V o )(0, y, 8). 

Since the leading term a(y) belongs to kerL(i(w, k)), and uj ^ 0, we have Tl s a — 0. Then the initial 
datum for N is 0(1) in H x (Jg, H^ 2 ): 



s-2\ 

y > 



A/"(0, y, 8) = U sai (y, 8) + eU s a 2 (y, 9) - J(V , V o )(0, y, 8) G H l (T e ,H, : 
This gives well-posedness over diffractive times: 

Proposition 3.2. The Cauchy problem (|3.20[) - (|3.21[) admits a unique solution on a time interval 
[0,T 2 *[ 7 with T 2 * > independent of e. Moreover, for all T < T 2 , we have the estimates : 

sup||(Vo,A0(r)|| H1(T H s-2 ) <C, 

[0,T] 

where the constant C = C(s,T) is independent of e and depends on s through the sum of norms 

Proof. The system (|3.20p being symmetric hyperbolic, local-in-time well-posedness is classical. 
Here the bounds are uniform in e. Indeed, the initial datum (|3.2ip is uniformly bounded with 
respect to e in Sobolev spaces i? 1 (Tg, H y ~ 2 ), s > 2 + 1/2. The L 2 estimate is uniform in e in 
spite of the large 1/e prefactor because the operator — ujdg + A(ei)kdg + L is skew-adjoint, and 
the commutator estimates are trivial because the operator —uidg + A(e\)kdg + Lq has constant 
coefficients. □ 

Back to the variable (Vo, N), where Vo and N are introduced in (|3.4j) and (|3.16[) respectively, by the 
definition of (Vo, A/"), we have that (Vo, N) is well defined over the long time interval [0, T 2 /e[. Back 



21 



to the variable V introduced in (|1.10[) . this gives well-posedness over the time interval [0,T 2 */e[ 
with the estimates for any T < T 2 * : 

sup\\n V(t)\\ HH s- 2) <C, sup \\Il s V(t)\\ HH m-*) < Ce, (3.22) 

[0,T/e] ^ [0,T/e] k ' V ' 

for the same constant C = C(s,T) as in Proposition 13.21 

We proved the long time existence of order 0(1/ e) of the solution of the Cauchy problem (|1.20[) . 
Back to the original variable v(t,y) = V(t, j/, 9) L_ -wt+ky , the following corollary follows immedi- 
ately. 

Corollary 3.3. Over the time interval [0,T 2 */e[ ; the Cauchy problem (|1.7p - (|1.9[) admits a unique 
solutionv oftheformv(t,x) = V(t,x, ky - — ), withV(t,x,9) G L°°([0, T/e] t , H 1 ^, H s - 2 CMl))) 

£ y 

for any T < T 2 *. Moreover, by Sobolev embedding, we have < C, where C = C(s,T) is given 

in Provosition \3.2[ 

In the next section, we consider the stability of the WKB solution from Section 2 and show two 
convergence results. 



4 Error estimates 

For the WKB solution V a from Section 2, we use the projection, rescalling and normal form 
reduction from Section 3: 

(V a ,V s a ) = (IL V a ,IL s V a ), W* = ^f, N a =W?-J(V a ,V a ), (VS,Af a )(r) = (V a ,N a )(r/e). 
where J is defined by (f37TT|) and (|3~18| . Then the system in (Vo,7V a ) is 

d T Vg + ^A V£ - ^-V a = 21T B(V Q , N a + J(V Q , V a )) + en i?, 

" B £ (4-1) 

d - N<1 ~ -T^ a = -2J(n B(V a , N a + J(V a , V a )), V a ) - 2eJ(V a , U R) + U S R 

As shown in Section 2.2, an existence time for (|4.1[) is T* (introduced in Section 2 as an existence 
time for (|2.41[1 ). Define the perturbations 

<P(T,y,0) = (V o -VS)(T,y,e) 1 ^(T,y,0) = (N-U a )(r,y,9), (4.2) 

then the couple ($, 'J) solves the following system over time interval [0, T£ 2 \ with T£ 2 = minlT-j*, T 2 } 
(T 2 * is introduced in Proposition 13 . 21 as an existence time of (Vo,AT) for (|3.20[) - (|3.21[1 ): 

cud 

d T <P + - = 2n B(V a , ¥) + i/ (V Q , N a , $, *) - en R, 

6 £ (4.3) 

- = H s (Vq,M, *) + 2eJ(V Q , U R) - IL S R, 

where Hq and H s are defined by 

Hq(Vq ,Af a , $, ¥) := 2U B(^,JV a + * + J(V a + $, V a + *)) + 2n B(V a , J($, $ + 2V a )), 
^s(Vq , A/" a , $, *) : =2 J(n B(V Q , AT + J(V a , V a )), V a ) 

- 2j(n s(v a + *, jv a + * + j(v a + v a + $)), v a + 

(4.4) 
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By (|2.46p . the initial datum for ($, is 

$(0) - (Vb - ^ Q )(0) - -en 6 - e 2 n 6!, 

*(o) = -(Vi - v; o )(0) - J(Vb, v )(o) + J(^ a , F a )(o) = -n s 6 + eb 2 , N " > ' 

where 6 and 6 X are defined by (|2~4Tj) and gUJ), and b 2 £ i? 1 ^, #*~ 2 ) is defined by 

62 = -II.61 + J(b + eh,Va(0)) + J(V(0),b + efei). 
We give two estimates for ($, ^) in the following two sections. 

4.1 First error estimate 

Here we assume only the polarization condition a(y) £ ker L(i(cj, k)). Then by (I4.9[) . we have 
$(0) = 0(e), #(0) = O(l) in H l (T e ,H^- 2 ). Then for the symmetric hyperbolic system flU), we 
have the following proposition: 

Proposition 4.1. T/ie Cauchy problem (|4.3I) - (|4.9I) admits a unique solution on time interval [0, T3 [ 
wt/j T 3 * > Ti 2 , and the following estimates hold for any T < T 3 : 

ll $ lli~([0,T],Hi(T9,Hr 2 )) ^ C '( £ + T )' ll*llL~([0,T],iyi(T9,Hj- 2 )) ~ C ' 

where the constant C = C(s,T) is independent of e and depends on s through the sum of norms 

M\h° + IK|| H i (Xj , ;H j-i 3 + IMlH^Hi-'y 

Proof. As in the proof of Proposition [321 the data being bounded in e, we have existence, unique- 
ness, and uniform bounds in e for short times with existence time T 3 independent of e. 
The source terms of (|4. 3[) is O(l), then we have the estimate: 

By the first equation of (|4.9p . the initial datum $(0) = 0(e). The classical H s estimate then gives 



Since we already know that ($,*) defined by (|4.2j) solves (|4.3l) - (14.9j) over time interval [0,7^ [ 
where T* 2 := min{Ti ,T 2 *}, we have T 3 * > Tf 2 by uniqueness. □ 

Remark 4.2. _ffere we s/ioiu first that the exact solution exists over long times 0(1/ e) before 
showing that it is approximated by the WKB solution on the intersection of their intervals of 
existence. Indeed, the Maxwell-Bloch structure allows us to find the normal form of the nonlinear 
equations, in contrast to, e.g., the linear normal form reduction o/[5]. 

Remark 4.3. For this first error estimate, we do not need the special choice of f given by (I2.43[) . 
The estimate (|4.6[) , as well as the first estimate (| 1 . 24[) in Theorem \l.3l hold for any regular function 
f that satisfies the transport equation (|2.38j) . e.g. f = 0. (Recall, f is introduced in Section 2.1 as 
a building block of the first corrector V 1 in WKB approximation. See (|2.18|) ). 

Back to the original time and variables, we immediately obtain the following corollary: 

Corollary 4.4. For any T < T^, we have the error estimate 

r \\Mv -v a )\\ L „ {[0>T/E]<HWi - 2)) <c(e + n 

\ ||n a (V - V a )\\ L ^([v,T/E],m(T e ,H s y - 2 )) ^ Ce - 



The constant C = C(s,T) is as in Proposition 14.1 
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4.2 Second error estimate 

In this section, we show the stability of the WKB solution in the following sense: for some T > 
independent of e, we have 

\\Mv niL~ ([ o,T/ e]TO ,^- 2 )) = °( £ )> W U '(V - vnh^ao.T/eiHHT^Hr 2 )) = °( £2 )' 

provided the initial perturbation 

\\U (V - V a )(0)\\ H1(TgMr i } = 0(e), \\U S (V - V a ){0)\\m<?e,H}-*) = °^)- ( 4 - 10 ) 

Here we assume that the initial corrector a\ satisfies (I2.51|) . Then by Lemma (12.11) . 

||n a (y - V a )\\ L ^ {[0T/£lH1{i:f)H s-2 )) = 0(e 2 ), 

and the initial datum of ($, (introduced in (|4.2p ) is 

$(0) = -£n 6-e 2 n 6i = 0(e), *(0) = -IL 8 b + eb 2 = eb 2 = 0(e). (4.11) 

For the remainder R (introduced in (|2.46[) and explicitly given in (|2.49p ). with the choice of / 
in (J213), by Lemma O we have that U S R = ~eR s with R s = 2B{V X ,V 2 ) + sB(V 2 , V 2 ) G 
L°°([0,T/e],H 1 (Te,H*- 2 )) for any T < Tj*. We now consider the rescaled variables 

$i = $/e, *i = (4.12) 

then the equation in ($i, <£~i) is 

a r $! + -1^0*1 - = 2n 0J B(v Q ,*i) + -#o(v ,^v$i,^i)-n OJ R, 

£ 5 " £ (4-13) 

9 r *i - = -iJ s (V a ,A^,e$i,e*i) +2J(V Q ,n i?) - fl s , 

with the initial datum 

$i(0) = -Hoft-eDoftx, *i(0) = 6 2 , (4.14) 

and where (Hq,H s ) are defined in (|4.4|l . Note that by bilinearity of £?, and pointwise bounds for 
the approximate solution, we have uniform bounds 

i|#o(V Q ,A/' a , £$,£*)) < C(|$| + |*|), J|i7 s (V a ,AA a , e$,e*)| < C(|$| + |*|). 

It is now classical to deduce uniform bounds for $i and *i in times 0(1): 

Proposition 4.5. With the choice of a x in (I2.51[) , £/ie Cauchy problem (|4.13[) - (I4.14[) admits a 
unique solution ($!,*!) on [OjTJf loii/i T4 > independent of e, and for any T < T£ : 

ll $ l|li~([0,T],fl'l(T e ,HJ- 2 )) - C ' H*lllL~([0,T],_H-i(T 8 ^r 2 )) ~ ^' ( 4 - 15 ) 

where the constant C = C(s,T) is independent of £ and depends on s through the sum of norms 
IHIff^ + ll a illiri(T 9 ,.fr»- 1 ) + ll a 2llHi( Te ^=- 2 )- 

Remark 4.6. For i/iis second error estimate, we do need the special choice of f in (|2.43[) . Together 
with the choice of initial corrector a\ in (|2.51[) . it allows us to rescale the solution (|4.12[) to obtain 
the desired estimates. 
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Back to the original time and variables, we immediately deduce the following corollary: 
Corollary 4.7. Let T* — min-fT^, T£ }, then for any T < T* , we have the error estimates 

WW -V a )\\ Loo{l ^ T/ElH1{Te ^ y - 2)) <Ce, \\Ils(V-Vn\\ L ~ il0 , T/elHHTe , H s-z }) <Ce^ (4.16) 
where the constant C = C(s,T) is as in Provosition \4:.5\ 

This shows stability, that is, the WKB approximate profile stays close to the exact profile over its 
existence time with an error estimate that is comparable to the initial error. 

4.3 Proof of Theorem [L3] 

We now sum up and prove Theorem II .31 

First, by Corollarv l3.3[ the function v{t, y) — V(t, y, 6)\ e= k V -^t solves (|1.7[l - (|1.9[) over time interval 
[0,T*[, and this gives the first result of Theorem ll.3l 

Second, by the Corollary 03] and Sovolev embedding if^Te, #*" 2 ) c L°°(T x R y ), the error 
estimates (| 1 . 24f) and (|1.26|) of Theorem 11.31 follow immediately. 

Finally, by the Coroll ary and Sovolev embedding if 1 (TV, if*" 2 ) C L°°{T g x R y ), we have the 
second error estimate f|l .25[) . This shows stability of the approximate solution. 
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